Exclusive electroproduction of <p mesons at 4.2 GeV 



K. Lukashin, 1 - 2 * E.S. Smith, 2 G.S. Adams, 28 E. Anciant, 7 M. Anghinolfi, 15 B. Asavapibhop, 20 G. Audit, 7 
T. Auger, 7 H. Avakian, 14 J. Ball, 3 S. Barrow, 13 M. Battaglieri, 15 K. Beard, 18 M. Bektasoglu, 25 W. Bertozzi, 21 
N. Bianchi, 14 A. Biselli, 28 S. Boiarinov, 17 B.E. Bonner, 29 S. Bouchigny, 2 D. Branford, 11 W.J. Briscoe, 16 
W.K. Brooks, 2 V.D. Burkert, 2 J.R. Calarco, 22 D.S. Carman, 24 B. Carnahan, 6 L. Ciciani, 25 R. Clark, 5 P.L. Cole, 32 ' 2 
A. Coleman, 34 t D. Cords, 2 P. Corvisiero, 15 D. Crabb, 33 H. Crannell, 6 J. Cummings, 28 P.V. Degtiarenko, 2 
L.C. Dennis, 13 E. De Sanctis, 14 R. DeVita, 15 K.S. Dhuga, 16 C. Djalali, 31 G.E. Dodge, 25 J. Domingo, 2 
D. Doughty, 8 - 2 P. Dragovitsch, 13 S. Dytman, 27 M. Eckhause, 34 H. Egiyan, 34 K.S. Egiyan, 35 L. Elouadrhiri, 8 ' 2 
A. Empl, 28 R. Fatemi, 33 R.J. Feuerbach, 5 J. Ficenec, 1 K. Fissum, 21 T.A. Forest, 25 A. Freyberger, 2 H. Funsten, 34 

S. Gaff, 10 M. Gai, 9 G. Gavalian, 35 * S. Gilad, 21 G. Gilfoyle, 30 K. Giovanetti, 18 P. Girard, 31 K.A. Griffioen, 34 
M. Guidal, 26 M. Guillo, 31 V. Gyurjyan, 2 J. Hardie, 8 ' 2 D. Heddle, 8 ' 2 F.W. Hersman, 22 K. Hicks, 24 R.S. Hicks, 20 
M. Holtrop, 22 J. Hu, 28 C.E. Hyde- Wright, 25 M.M. Ito, 2 D. Jenkins, 1 K. Joo, 33 § J. Kelley, 10 M. Khandaker, 23 
K. Kim, 19 K.Y. Kim, 27 W. Kim, 19 A. Klein, 25 F.J. Klein, 2 ** M. Klusman, 28 M. Kossov, 17 L.H. Kramer, 12 ' 2 
5. Kuhn 25 J.M. Laeet. 7 D. Lawrence 20 A. Lonehi 6 J.J. Manak 2 tf C. Marchand. 7 S. McAlecr 13 J. McCarthy 



M. Holtrop, 22 J. Hu, 28 C.E. Hyde- Wright, 25 M.M. Ito, 2 D. Jenkins, 1 K. Joo, 33 § J. Kelley, 10 M. Khandaker, 
K. Kim, 19 K.Y. Kim, 27 W. Kim, 19 A. Klein, 25 F.J. Klein, 2 ** M. Klusman, 28 M. Kossov, 17 L.H. Kramer, 12 ' 
S.E. Kuhn, 25 J.M. Laget, 7 D. Lawrence, 20 A. Longhi, 6 J.J. Manak, 2 tt C. Marchand, 7 S. McAleer, 13 J. McCarthy, 
J.W.C. McNabb, 5 B.A. Mecking, 2 M.D. Mestayer, 2 C.A. Meyer, 5 K. Mikhailov, 17 R. Minehart, 33 M. Mirazita, 14 
R. Miskimen, 20 V. Muccifora, 14 J. Mueller, 27 G.S. Mutchler, 29 J. Napolitano, 28 S. Nelson, 10 G. Niculescu, 24 
I. Niculescu, 16 B.B. Niczyporuk, 2 R.A. Niyazov, 25 A. Opper, 24 G. O'Rielly, 16 J.T. O'Brien, 6 K. Park, 19 
K. Paschke, 5 E. Pasyuk, 3 G.A. Peterson, 20 S. Philips, 16 N. Pivnyuk, 17 D. Pocanic, 33 O. Pogorelko, 17 E. Polli, 14 

S. Pozdniakov, 17 B.M. Preedom, 31 JW. Price, 4 L.M. Qin, 25 B.A. Raue, 12 ' 2 A.R. Reolon, 14 G. Riccardi, 13 
G. Ricco, 15 M. Ripani, 15 B.C. Ritchie, 3 F. Ronchetti, 14 P. Rossi, 14 D. Rowntrec, 21 P.D. Rubin, 30 F. Sabatie, 25 
K. Sabourov, 10 CW. Salgado, 23 ' 2 V. Sapunenko, 15 R.A. Schumacher, 5 V. Serov, 17 Y.G. Sharabian, 35 § J. Shaw, 20 
S. Simionatto, 16 A. Skabelin, 21 L.C. Smith, 33 D.I. Sober, 6 A. Stavinsky, 17 S. Stepanyan, 35 « P. Stoler, 28 
I.I. Strakovsky, 16 M. Tahiti, 15 S. Taylor, 29 D. Tedeschi, 31 ' 2 R. Thompson, 27 M.F. Vineyard, 30 A. Vlassov, 17 
K. Wang, 33 H. Weller, 10 L.B. Weinstein, 25 R. Welsh, 34 D.P. Weygand, 2 S. Whisnant, 31 E. Wolin, 2 L. Yanik, 16 
A. Yegneswaran, 2 J. Yun, 25 J. Zhao, 21 B. Zhang, 21 Z. Zhou 21 « 

(The CLAS Collaboration) 

1 Virginia Polytechnic Institute and State University, Department of Physics, Blacksburg, VA 24-061, USA 
2 Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA 
3 Arizona State University, Department of Physics and Astronomy, Tempe, AZ 85287, USA 
4 University of California at Los Angeles, Department of Physics and Astronomy, Los Angeles, CA 90095, USA 
5 Carnegie Mellon University, Department of Physics, Pittsburgh, PA 15213, USA 
6 Catholic University of America, Department of Physics, Washington, DC 20064, USA 
7 CEA Saclay, DAPNIA-SPhN, F91191 Gif-sur-Yvette Cedex, France 
8 Christopher Newport University, Newport News, VA 23606, USA 
9 University of Connecticut, Physics Department, Storrs, CT 06269, USA 
10 Duke University, Physics Department, Durham, NC 27706, USA 
11 Edinburgh University, Department of Physics and Astronomy, Edinburgh EH9 3JZ, United Kingdom 
12 Florida International University, Department of Physics, Miami, FL 33199, USA 
13 Florida State University, Department of Physics, Tallahassee, FL 32306, USA 
14 Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali di Frascati, CP. 13, 00044 Frascati, Italy 
15 Istituto Nazionale di Fisica Nucleare, Sezione di Genova e Dipartimento di Fisica dell'Universita, 16146 Genova, Italy 
The George Washington University, Department of Physics, Washington, DC 20052, USA 
17 Institute of Theoretical and Experimental Physics, Moscow, 117259, Russia 
18 James Madison University, Department of Physics, Harrisonburg, VA 22807, USA 
19 Kyungpook National University, Department of Physics, Taegu 702-701, South Korea 
20 University of Massachusetts, Department of Physics, Amherst, MA 01003, USA 

21 M. I. T. -Bates Linear Accelerator, Middleton, MA 01949, USA 
22 University of New Hampshire, Department of Physics, Durham, NH 03824, USA 



24 



23 Norfolk State University, Norfolk, VA 23504, USA 
Ohio University, Department of Physics, Athens, OH 45701, USA 



2s 



25 Old Dominion University, Department of Physics, Norfolk, VA 23529, USA 
26 Institut de Physique Nucleaire d'Orsay, IN2P3, BP 1, 91406 Orsay, France 
27 University of Pittsburgh, Department of Physics, Pittsburgh, PA 15260, USA 
Rensselaer Polytechnic Institute, Department of Physics, Troy, NY 12181, USA 



29 Rice University, T.W. Bonner Nuclear Laboratory, Houston, TX 77005-1892, USA 
30 University of Richmond, Department of Physics, Richmond, VA 23173, USA 
31 University of South Carolina, Department of Physics, Columbia, SC 29208, USA 



1 



University of Texas, Department of Physics, El Paso, TX 79968, USA 
33 University of Virginia, Department of Physics, Charlottesville, VA 22903, USA 
College of William and Mary, Department of Physics, Williamsburg, VA 23185, USA 
35 Yerevan Physics Institute, 375036 Yerevan, Armenia 
(February 7, 2008) 



Abstract 

We studied the exclusive reaction ep — ► e'p'<f) using the <f> — > K + K~ decay mode. The data were 
collected using a 4.2 GeV incident electron beam and the CLAS detector at Jefferson Lab. Our 
experiment covers the range in Q 2 from 0.7 to 2.2 GeV 2 , and W from 2.0 to 2.6 GeV. Taken 
together with all previous data, we find a consistent picture of cj> production on the proton. Our 
measurement shows the expected decrease of the f-slope with the vector meson formation time cAt 
below 2 fm. At (cAr) = 0.6 fm, we measure = 2.27 ±0.42 GeV -2 . The cross section dependence 
on W as W° 01 at Q 2 = 1.3 GeV 2 was determined by comparison with <j> production at HERA 
after correcting for threshold effects. This is the same dependence as observed in photoproduction. 

PACS number(s): 13.60.Le, 13.60.-r, 12.40.Vv, 14.40.Cs 
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I. INTRODUCTION 

Vector meson photo- and electroproduction have been 
important tools used to understand the hadronic prop- 
erties of the photon For low values of the four- 
momentum transfer squared, the photon interacts with 
the target predominantly through vector meson interme- 
diate states that diffractively scatter off the target. This 
process, shown in Fig.|l|a, was originally described within 
the framework of the Vector-Meson Dominance (VMD) 
model. The production of the <j) meson through this 
mechanism may be interpreted in terms of the hadronic 
structure of the photon, that couples to a virtual meson 
with a strength proportional to the square of the charge 
of its constituent quarks. Due to the dominant ss quark 
component in the <f> meson, quark exchange (e.g. meson- 
exchange) mechanisms and s-channel resonance produc- 
tion are strongly suppressed (|-||]. As a consequence, 
4>p scattering at low four-momentum transfer proceeds 
primarily through pomeron exchange, similar to hadron- 
hadron diffractive scattering. 



b) 




FIG. 1. Representation of (f> production by: a) the 
VMD model, and b) the Donnachie and LandshofF 
pomeron-exchange model. 

It is generally believed that the underlying mechanism 
for pomeron exchange is multi-gluon exchange, where 
the simplest possibility requires at least two gluons since 
all hadrons are color singlets. A simplification to these 
calculations was introduced by Donnachie and Landshoff 
HI?), who proposed a model whereby the pomeron cou- 
ples to quarks inside the interacting hadrons as shown 
in Fig. [l]b. Calculations within this context have been 
applied to <j> electroproduction to study the quark sub- 
Mi and to photoproduction at large 
In these models the cross 
section increases slowly with center-of-mass energy, IV, 
reflecting the pomeron trajectory. 

At high Q 2 the pomeron can be resolved into two- 
gluon exchange, and predictions for hard diffractive elec- 
troproduction of vector mesons can be made within the 
context of perturbative QCD (O] . At lower energies (IV 



structure of mesons 
momentum transfer 



~ 10 GeV), quark exchange mechanisms [[l3],fL4) become 
significant for the production of vector mesons with va- 
lence u and d quarks, but play a limited role in the pro- 
duction of 4> mesons. 

The hadronic structure of the photon arises from fluc- 
tuations of the virtual photon into short-lived quark- 
antiquark (qq) states of mass My during a formation 
time 01 



At = 



2v 



(Q 2 



Ml)' 



(1) 



where -Q 2 is the squared mass and v is the laboratory- 
frame energy of the virtual photon (see Appendix A for 
notation) . The effect of the formation time on the prop- 
agation of these virtual quantum states in strongly inter- 
acting media has been observed for p mesons propagat- 
ing inside a proton |l5|] and inside nuclear targets [ ^6[ . 
To date, no clear dependence on the formation time has 
been observed in <j) meson production by virtual photons 

This paper presents measurements of exclusive <fi me- 
son electroproduction off a proton target for 2.0 < 
W < 2.6 GeV and 0.7 < Q 2 < 2.2 GeV 2 where there 
is extremely limited data. In this kinematic regime, 
the short formation distance ^ of the virtual qq state 
(0.35 < cAt < 0.75 fm) limits the time for interaction 
and probes the <p production mechanism at small forma- 
tion times. 

In Section H we present the details of our experimental 
techniques and data analysis. It concludes with values 
for the measured t-slopes and total cross sections. In 
Section [II we compare our results with previous data, 



and compare with a geometrical model of the relation 
between formation time and t-slope. The model is dis- 
cussed in some detail in Appendix |b[ 



II. EXPERIMENT 

The experiment was performed using the CEBAF 
Large Acceptance Spectrometer (CLAS) [|l9],^0) in Hall 
B of the Thomas Jefferson National Accelerator Facility 
(Jefferson Lab). The data were taken with a 4.2 GeV 
electron beam incident on a 5.0 cm liquid hydrogen tar- 
get in March and April of 1999. The CLAS torus magnet 
current was set to 2250 A, bending negatively charged 
particles toward the beam axis. The trigger required 



x In the literature the formation distance is also referred to 
as coherence length. 
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a single scattered electron signal, identified as a coinci- 
dence of the forward electromagnetic calorimeter (EC) 
PH and Cerenkov counters [^]. Data were recorded at 
an instantaneous luminosity of 0.6 • 10 34 cm~ 2 s _1 and 
a typical livetime of 95%. This data set has a live-time 
corrected integrated luminosity Li nt — 1.49 • 10 39 cm~ 2 . 

Data Reduction. In order to reduce the data sample 
to a manageable size, the data of interest were first pre- 
selected using very loose requirements on particle iden- 
tification, missing mass, and the requirement for W to 
be above 1.8 GeV. The tj> mesons were identified through 
their K + K~ decay mode. Because of the small accep- 
tance of K~ due to the CLAS magnetic field setting, we 
required only three final-state particles to be detected: 
electron, proton, and K + . The K~ was reconstructed 
by identification in the epK + (X) missing mass. The mo- 
menta of charged tracks were reconstructed from their 
curvature in the CLAS magnetic field using a system of 
drift chambers ||| . The data reduction process selected 
about 82,000 events for further analysis. The size of this 
filtered data sample was compact (ss 0.5 GigaByte) and 
easily manageable in comparison with the size of the en- 
tire data set (« 1 TeraByte). 




FIG. 2. a) Electron momentum versus total deposited en- 
ergy in EC. The solid lines show the applied cuts, b) Energy 
deposited by the TOF-identified 7r~'s in the outer EC layers 
versus energy deposited in the inner EC layers. The solid 
line shows the applied cut Ei n > 0.04 GeV, which retains all 
good electron candidates. 

Electron Identification. In addition to a fiducial 
requirement that an electron hit be at least 10 cm from 
the outer edge of the electromagnetic calorimeter, cuts 
on energy deposition in the EC were applied in order to 
avoid misidentification of tt~ as e~ . The total energy 
deposited by an electron in the EC is proportional to 
the momentum determined by magnetic analysis. This 
dependence is illustrated in Fig. ||a. The electron band 
with the width of the EC resolution is clearly seen. In 
order to cut out hadronic background, we applied cuts 
around this band (the solid lines in Fig. 0a) . An addi- 



tional improvement in e~ /ir~ separation was achieved 
by cutting out the ir~ signal based on the energy de- 
posited in the inner layer of the calorimeter as shown in 
Fig. ||b. The cluster of entries to the left of the line is 
the ir~ signal in the EC. The solid line is the applied cut 
(Ei n > 0.04 GeV) to eliminate pions. To determine this 
cut we used n~ identified by the Time-of-Flight system 
(TOF) of the CLAS H). 

Hadron Identification. The identification of 
charged hadrons is illustrated in Fig. ||. The distribu- 
tion of positively charged particle momenta versus re- 
constructed mass is shown in Fig. ||a. Proton, kaon, 
and positive pion bands are clearly distinguished. The 
width of the reconstructed mass increases with momen- 
tum. However, there is no systematic dependence after 
careful timing calibration of the detector p4]-p6[ . 

K + Identification. In order to optimize the signal- 
to-background ratio in kaon identification, the kaon mo- 
mentum range was divided into six bins. In each bin 
the mass distribution was fitted to a Gaussian with a 
polynomial background to determine the characteristics 
of the K + peak. An example of this procedure is shown 
in Figs. ||a-b. The horizontal lines in Fig. [|]a show the 
momentum bins for K + identification, and the fitting 
result for one of the bins is illustrated in Fig. |3)h> . To 
identify kaons, ±2cr cuts were applied around the mean 
value (m K +). 




0.5 1 1.5 0.3 0.4 0.5 0.6 0.7 



M T0F (GeV) M TOF (GeV) 

FIG. 3. a) Positively charged particle momentum versus 
reconstructed mass for the pre-selected event sample. The 
apparent separation between kaons, pions and protons at high 
momenta is due to the data pre-selection cuts. The horizon- 
tal lines show the binning in kaon momenta; b) K + recon- 
structed mass distribution in the momentum bin from 0.9 to 
1.2 GeV. The background is due to pion misidentification. 

Proton Identification. The proton signal is very 
clean and does not have a significant background con- 
tribution. For proton identification we applied a simple 
reconstructed mass cut from 0.8 to 1.1 GeV. 

K~ Identification. We identified the K~ using the 
missing mass technique. The K~ band is clearly seen 
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in Fig. [|a. The selection used ±2<r cuts around the K~ 
peak. The invariant mass, Mkk, of the K + K~ is com- 
puted using the known mass of the kaons, the measured 
momentum of the K + and the missing momentum of 
the event for the K~ . We note that because the masses 
are large compared to the momenta of the particles, this 
quantity has significantly better resolution than the epX 
missing mass. 

Identification of the Signal. Applying the electron 
and hadron identification cuts described above, we iden- 
tified about 3800 events of the epK + K~ final state. In 
order to eliminate events caused by false triggers on low 
energy electrons (e.g. from ir° Dalitz decays) we also 
required the energy transfer, v — E e — E e i , to be smaller 
than 3.5 GeV. The selected sample includes 4> mesons, 
high mass hyperons, and background events that come 
from particle misidentification. 
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FIG. 4. The (f> channel separation technique, a) epK + X 
missing mass versus epX missing mass. The horizontal 
lines show the selection of K~ . b) Mkk mass spectrum of 
epK + K~ events, c) eK + X missing mass distribution. The 
line shows the A(1520) cut. d) Mkk mass distribution with 
the A(1520) cut applied. 

The most important features of the final selection are 
shown in Figs.[|a-c. In the scatter plot of epK + X versus 
epX missing mass (see Fig.|^a) the signal of the epK + K~ 
final state is clearly distinguished from the rest of the 
data. The lines solid show the ±2cr selection cuts in the 
reconstructed K~ mass. Fig. ||b shows the Mkk mass 
distribution of the selected final state with a prominent 



peak due to exciting <f> particles. To extract the total <p 
yield, we fitted the peak with a Gaussian (the integral is 
shown as the filled area in the plot) and the background 
with an empirical phase space function, 



f(M KK ) = Ay M KK - Ml + B(M KK - Mf h ) , (2) 

where the threshold M th = 0.987 GeV. The fit gives 
= 248, a mean value (M KK ) = 1019.1 ± 0.6 MeV, 
and a = 6.0 ± 0.6 MeV, where the width of the peak is 
dominated by the resolution of CLAS. ^ The </> signal-to- 
background ratio is 0.7 within ±2er from the mean value 
of the 4> peak. 

The primary source of physical background consists 
of high mass hyperons, ep — > e'K + Y*, with a subse- 
quent decay Y* — > NK. The production and decay 
amplitudes of these particles are not well known. The 
main channel is the A(1520) with a cross section larger 
than <f> production. Additional contributions come from 
A(1600), A(1800), A(1820), £(1660), and £(1750), which 
have large branching ratios for decay into the NK chan- 
nel (27) . These backgrounds were investigated by Monte 
Carlo using exactly the same algorithms as the experi- 
mental data in order to optimize selection cuts. In order 
to minimize the number of A(1520) in the data sample, 
we require Mx(eK + X) to be greater than 1.56 GeV. 
The cut is shown for the data sample with the solid line 
in Fig. He. 
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FIG. 5. Kinematic distributions of the selected (f> events: 
a) Q 2 versus W. b) Ar versus Q 2 . The dashed lines indicate 
the binning, used later, in Q 2 and Ar; the solid lines show 
the range of W used in the analysis. 

The Mkk mass distribution with the A(1520) cut ap- 
plied is shown in Fig.^d. The simultaneous fit of the cf> 
peak and the background gives — 197, a mean value 



2 The mass of the (j> is 1019.417 ± 0.014 MeV, and the decay 
width (FWHM) is 4.458 ± 0.032 MeV @- 
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(M KK ) = 1019.4±0.9 MeV, and a = 6.4±1.1 MeV. The 
4> signal-to-background ratio is improved, and equals 1.3 
within ±2(7 of the 4> peak. The remaining background, 
consistent with phase space, is due to high-mass hyperon 
states, non-resonant K + K~ production and experimen- 
tal misidentification of a ir + as a K + (events under the 
K + peak in Fig.^b). We note that the level of the back- 
ground under the (f> peak depends on the fitting proce- 
dure, and will be addressed when we discuss systematic 
errors. 

The kinematic range of the data sample is shown in 
Fig. |. The range of Q 2 varies from 0.7 to 2.2 GeV 2 , 
W from 2.0 to 2.6 GeV, and At from 1.8 to 4.0 GeV" 1 
(cAr from 0.35 to 0.79 fm). The small values of cAt 
indicate that the formation distance in our kinematic 
regime is below the hadron size, 2 w 2 fm. The data 
binning to calculate the exponential i-slope (see below) 
is indicated in Fig. pi by horizontal dashed lines, which 
show the ranges of Q 2 (integrated over At) and At (in- 
tegrated over Q 2 ). In both cases the data range in W 
is the same (solid lines in Fig. [s|a). We note that finer 
binning in Q 2 and W is used for the evaluation of the 
cross sections integrated over t'. 




(GeV) 
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FIG. 6. Side-band background subtraction technique, a) 
Location of the side bands; b), c) and d) distributions of 
events in the signal region (histograms) and in the side-bands 
(open squares) versus Q 2 , W and | t — t m in • 

Ideally, with enough statistics and an understanding 
of the background shape, fits would be used to extract 
the signal yield in every kinematic bin of interest. With 



limited statistics this is not possible, and we proceeded 
by using a side-band subtraction technique. 

Background Subtraction. The side-band tech- 
nique, as illustrated in Fig. ^, was used to determine 
the background distribution as a function of Q 2 , W and 
—if . The signal region was determined within a ±2cr cut 
around (M^) after excluding the A(1520) from the final 
state data sample. The side bands were located ±3<r 
away from the (j) peak, and the number within the band 
was scaled to the background as determined by the fit 
(see Fig.||d). The normalized side-band events were then 
subtracted in each distribution of interest. This proce- 
dure is illustrated in Fig. ^ for the entire data set and 
was repeated for each kinematic region defined in Table 
1. 
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region 


Q 2 and cAt 
range 


(Q' 2 ) 

(GeV 2 ) 


(cAr) 
(fm) 


b,p 

(GeV" 2 ) 


All Data 


0.7 - 2.2 GeV^ 
0.35 - 0.75 fm 


1.02 


0.6 


2.27±0.42 


Low Q 2 


0.7 - 1.2 GeV^ 


0.87 




2.31±0.59 


High Q' 2 


1.2 - 2.2 GeV 2 


1.47 




2.10±0.52 


Low cAr 


0.35 - 0.55 fm 




0.49 


2.04±0.42 


High cAr 


0.55 - 0.75 fm 




0.63 


2.12±0.46 



TABLE 1. The measured values of the t-slope parameter, 
&0, fitted to the data for —if < 1.2 GeV 2 . The errors are 
statistical only. 
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FIG. 7. Acceptance as a function of Q 2 and —t'. 

Acceptance. For the calculation of the acceptance, 
we used a GEANT-based simulation of CLAS, taking 
into account trigger efficiency, problematic hardware 
channels, and the CLAS resolution. The Monte Carlo 
event sample was generated assuming the VMD model 
for (f> electroproduction. Two iterations in the accep- 
tance calculation were made to adjust the VMD param- 
eters to be close to the data. In each kinematic region, 
the acceptance was calculated from the ratio of recon- 
structed to generated <p events with the same kinematics 
and particle identification cuts that were applied to the 
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data. Fig. fj] shows the acceptance as a function of Q 2 
and — t' for the entire data set. This procedure was also 
used to calculate the acceptance as a function of W and 
At in each kinematic bin. 




b) 



Q 2 (GeV*) 

FIG. 8. Simulated data: a) <j> peak, convoluted with the 
measured CLAS resolution, with radiative effects turned off 
(solid), and turned on (hatched histogram), b) Inverse radia- 
tive correction factor, 1/F rad , as a function of Q 2 . 

Radiative Corrections. For the calculation of the 
radiative corrections, we used the peaking approxima- 
tion p8[ . We define the radiative corrections in each 
bin of every kinematic variable as the ratio F rad — 
N n0 rad/N rad , where N rad and N norad are the generated 
<fi yields with radiative effects turned on and off, respec- 
tively. The model for the <p production cross section 
employed for the computation of acceptance was also 
used for the studies of radiative corrections. The ratios 
were calculated with the same kinematics and particle 
identification cuts that were applied to the data. The 
simulated <j) mass distributions with and without radia- 
tive effects are shown in Fig. ||a. The inverse radiative 
correction factor, 1/F rad , as a function of Q 2 is shown 
in Fig. ||b. The correction factors as a function of —t' 
in all four kinematic regions are of the order of 1.4 and 
uniform over the kinematics considered here. 



Q 2 bin 
(GeV 2 ) 


{W) 
(GeV) 


(e) 


T{Q\W) 
(10- 4 GeV- 3 ) 


<t(Q*,W) 

(nb) 


0.8 - 1.0 


2.37 


0.51 


1.50 ± 0.15 


27.6 ± 6.1 


1.0 - 1.2 


2.31 


0.50 


1.12 ± 0.10 


24.2 ± 5.4 


1.2 - 1.4 


2.28 


0.49 


0.879 ± 0.067 


23.0 ± 5.2 


1.4 - 1.6 


2.28 


0.44 


0.701 ± 0.050 


20.8 ± 5.7 


1.6 - 1.8 


2.25 


0.42 


0.562 ± 0.033 


14.5 ± 6.4 



TABLE 2. The averaged values of W, e, F(Q 2 ,W), and 
a(Q 2 ,W) as a function of Q 2 . The numbers given for the 
virtual photon flux, V(Q 2 ,W), computed event-by-event, are 
the mean and the standard deviation for the bin. 

Data Normalization. The final step in the analysis 
procedure was the normalization of the <f> yield to the 



integrated luminosity, the virtual photon flux, and all 
calculated corrections as: 



N^/BR^k+k- F acc F rad F wm 



(3) 



AQ 2 AW 27rT(Q 2 ,W)L mt 

where AQ 2 and AW are the bin widths in Q 2 and W, 
T(Q 2 ,W) is the virtual photon flux, L int is the inte- 
grated luminosity, N$ is the <j> yield in the bin, F acc is 
the acceptance factor in a given bin, F wm is a small 
correction factor for production from the target win- 
dows (w 1%), F rad is the radiative correction factor, 
and BR = 0.492 ± 0.007 is the decay branching ratio for 
</> — ► K + K~ ^7|. The virtual photon flux was calcu- 
lated on an event-by-event basis and averaged for each 
kinematic bin as 

W^^^XfS^-, (4) 



M P Q 2 



1 



where M p is the mass of the proton, E e is the electron 
beam energy, and e is the polarization of the virtual pho- 
ton: 



AE e {E e -v)-Q 2 
4E e (E e -v)+2v 2 + Q 2 



(5) 



Cross Section, <r(Q 2 ,W). The cross section inte- 
grated over all t' , a(Q 2 ,W), was extracted in five bins 
over a Q 2 range from 0.8 to 1.8 GeV 2 with a bin width 
of 0.2 GeV 2 . The range in W was determined as the 
allowed kinematic range for each Q 2 . The binning, val- 
ues of the virtual photon flux used during normalization, 
T(Q 2 ,W), and the measured cross section are given in 
Table 2. The table shows statistical errors only. 




2.5 3 

(GeV 2 ) 

FIG. 9. The da/dt' differential cross section for exclusive 
<j> electroproduction off the proton with fits to the entire —t' 
range (dashed) and — t' less than 1.2 GeV 2 (solid). 
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Differential Cross Section, da/dt'. The measured 
cross section, da/dt' ', is generally parameterized at small 
-t'by 



da 
dt 1 



A, 



J>d>t' 



(6) 



The entire t' range (0 < -t' < 2.6 GeV 2 ) can be fit to a 
single exponential with a slope b& = 1.61±0.31 GeV -2 
and a x 2 = 0.9/DF. However, Eq.||is only expected to be 
valid at small — t', so we have restricted our analysis to 
—t' less than 1.2 GeV 2 which also allows direct compari- 
son to previous measurements. For this restricted range, 
we obtain i> = 2.27±0.42 GeV" 2 (solid line in Fig.g). 
We also performed fits in the four overlapping kinematic 
regions specified in Table 1: two ranges in Q 2 (integrated 
over cAr), and two ranges in cAr (integrated over Q 2 ). 
The results of these fits are given in Table 1. 

We note that at larger —t', there is an apparent change 
in the slope of the distribution with a break at —t' « 1.3 
GeV 2 . This suggests that additional mechanisms may 
be present at — t' > 1 GeV 2 . Despite the fact that the 
break is not statistically significant, we discuss possi- 
ble mechanisms for a slope change. A similar pattern 
is observed in hadron-hadron elastic scattering [p9| , ^o[ , 
where a dip is observed at —t « 1.4 GeV 2 followed by a 
second maximum at —t « 1.8 GeV 2 . However, photo- 
production data do not show a change in the slope for 
-t < 2 GeV 2 0. S-channel production of resonances 
results in a large measured value of — t'. However, there 
are no known resonances that decay into <pN. Finally, 
imperfect background subtraction could also lead to an 
enhancement at large —t', but should be subsumed into 
our quoted systematic errors. 



Source 


Aa{%) 


A6 (%) 


Target stability 


0.7 




Target walls 


1.0 




Acceptance 


7.8 


5.0 


Radiative corrections 


4.7 




Background subtraction 


5.4 


4.6 


Total 


10.7 


6.8 



TABLE 3. Summary of the contributions to the systematic 
errors. 

Systematic Errors. Estimates of our systematic er- 
rors for the cross section, Act, and the r-slope parameter, 
A&0, are given in Table 3. The errors are averaged over 
the kinematics of the experiment, although the lowest Q 2 
cross section point may have about twice this systematic 
uncertainty due to the steepness of the acceptance func- 
tion (see Fig.^,). To estimate the systematic errors due 



to background subtraction, a complete analysis of the 
cross section and r-slope parameter was performed us- 
ing two different assumptions for the shape of the back- 
ground: phase space and a constant. The difference be- 
tween these results is quoted as the systematic error due 
to background subtraction. The systematic errors due 
to acceptance and radiative corrections are discussed in 
Ref. |3l| and |B2L respectively. Additional details can be 
found in Ref. p5[ . We note that the overall uncertainty 
is dominated by statistical errors. 



a 

6 



10 



10 



10 




- CLAS Lukashin, W » 2.3 GeV 

- HI Adloff, W » 75 GeV 

- ZEUS Derrick, W = 70, 94, 99" GeV 

- CORNELL Cassel, W » 2.7 GeV. 

- BONN Besch, W » 2.2 GeV (Q 2 =fy 
SLAC Ballam, W = 2.5, 4.3 GeV (Q* 



1 10 

(Q 2 + M/)(GeV 2 ) 

FIG. 10. The (f> meson cross section dependence on Q 2 for 
photo- and electroproduction. Electroproduction data HI 
Adloff are from Ref. [34], ZEUS Derrick from Ref. [35,36], 
and CORNELL Cassel from Ref. [15]. Photoproduction data 
BONN Besch are from Ref. [37], and SLAC Ballam from Ref. 
[38]. The solid and dashed curves are the pomeron-exchange 
model predictions for W=70 GeV and for 2.0 < W < 2.6 
GeV, respectively [10]. 



III. RESULTS 

Cross Section Dependence on Q 2 and W. The 

world data on elastic virtual photon production of <p 
mesons are shown as a function of Q 2 in Fig. [lO], and 
as a function of W in Fig. |ll|. Selected photoproduction 
data are also plotted for completeness. ^ We show the 
data on both plots with common symbols. 



' Additional data of (f> production on nuclear targets |33j] are 
available at (W) M 14 GeV. 
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All HERA data |34|-|36|] correspond to W ranging from 
40 to 130 GeV, where the gluonic density in the proton at 
low x = Q 2 /2M p is plays a significant role. Only the Cor- 
nell measurement [H| exists at low W, corresponding to 
x in the valence region. Q For the high-energy data, the 
Q 2 behavior of the cross section is well described by the 
vector meson propagator squared. The data are not yet 
in the asymptotic perturbative QCD regime where the 
longitudinal cross section for vector meson production 
is dominant, and should scale as Q~ 6 |3j|. Neverthe- 
less, the longitudinal contribution becomes increasingly 
important and must be treated systematically. For ex- 
ample, p mesons in muoproduction at large Q 2 are found 
to be dominantly in the helicity zero spin state E3] . 



- 1 



© 

D 
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10 



10 



; 0.0 


0.0 , 


o.o <>.--- " 

: 




- 0' 






2.7 ■ 

3.4 ■ 


13 W 


4.8 a 


r t/f 


7.5 m 

m 1 8.2 
12.1 ^ 

1 14.7 



10 



10 

W ( GeV ) 



FIG. 11. The 4> meson cross section dependence on W for 
photo- and electroproduction. The Q 2 values of the mea- 
surements are printed near the corresponding data points. 
All data points are from the same References as in Fig. 10. 
The curves, described in the text, correspond to Q 2 of 0, 1.3 
and 2.2 GeV 2 . 

Pomeron-exchange models, such as those described in 
the introduction, reproduce the Q~ 4 behavior of the data 
at large Q 2 . The predictions of a model |[(|, based on 
the Donnachie-Landshoff pomeron exchange (Fig. |l|b), 
are shown in Fig. [l^ for the W range of our experiment 
(2.0 < W < 2.6 GeV) and at W=70 GeV. The model 



4 We note that data points from Ref. [la] have different in- 
tegration ranges for the cross section as a function of Q 2 and 



describes the data reasonably well at high W, and re- 
produces the trend at low W, but overestimates the new 
cross section results presented here. We note that our 
data are close to the <j> production threshold, where the 
cross section increases rapidly as a function of center-of- 
mass energy. In the model of Pichowsky ^ , the transi- 
tion from a cross section that slowly decreases with Q 2 
to one that falls off as Q~ A occurs at a threshold that 
increases with the current-quark mass of the vector me- 
son. No clear threshold is visible in the <\> data, but the 
scarcity of points precludes drawing conclusions. 

The photoproduction cross section increases slowly 
with W, reflecting the pomeron trajectory. At higher 
Q 2 , a stronger dependence on W has been observed in 
preliminary analysis of HERA data ^dj . If the cross sec- 
tion is parameterized as W 5 , 8 varies from about 0.2 for 
photoproduction to 8 ~ 0.7 at a Q 2 of 8 GeV 2 . This in- 
creased dependence of the cross section on W has been 
interpreted as being due to the rise of the gluon momen- 
tum density in the proton at small x j3Sj ] . 

To be able to extract the W dependence by comparing 
our measurement at Q 2 — 1.3 GeV 2 to HERA data at 
the same Q 2 and (W) — 75 GeV, threshold effects must 
be taken into account. For example, threshold behavior 
can be clearly seen in the photoproduction data JO] (see 
Fig. |ll|). The reduced phase space near threshold be- 
haves as (Pcf>/q) 2 , where and q are the center-of-mass 
three-momenta of the <fi and virtual photon, respectively. 
This dependence of the cross section on W can be pa- 
rameterized as 



a(W) = G 



W 

Wo 



(7) 



W presented in Fig. ^ and Fig. [H]. 



Correcting for the threshold factor, our measurement of 
the cross section becomes <r C orr(Q 2 = 1.3) = 110±27 
nb, and using the HERA measurement, a(Q 2 — 1.3) = 
220±51 nb (34), we obtain 8 = 0.2±0.1. The quoted un- 
certainties were obtained by summing the statistical and 
systematic errors in quadrature. This slope is consistent 
with that measured in photoproduction. The curves of 
a(W) are shown in Fig.0 for Q 2 of 0, 1.3 and 2.2 GeV 2 , 
and 8 = 0.2. The curves are normalized to the HERA 
data ((To, W ) that are far from the production thresh- 
old. 

Dependence of t-slope on cAr. The dependence of 
the i-slope, 6^, on formation distance, cAt, for (j) meson 
production is shown in Fig. ^| together with previous 
data. In the terminology of the uncertainty principle, 
At is the time during which the virtual photon, with 
mass \/~Q 2 i can fluctuate into a <j> meson |Q . We expect 
that bt/) should decrease at low At as the interaction be- 
comes more point-like. The previous electroproduction 
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measurements |l^,|l7],|l8| do not show the expected be- 
havior. However, a consistent picture emerges when we 
include photoproduction data as well. Both of our data 
points (solid stars) lie in the region of cAt below 1 fm, 
and show a decrease of 6^ with decreasing formation time 
when combined with other data. This is consistent with 
the well-measured dependence for p meson production 
|fl5f as discussed in Appendix B. To fit the <f> meson data 
to Eq. B5, we constrain the parameter to the value 



extracted from the fit to the p data (Eq. B£ ) . This yields 
& (cAr) = (6.87 ± 0.17) [ 1 - e - cA -/ 2 (0-78) 1 (g) 



with x 2 /DF = 4.8. The fit to the (j> data is shown in 
Fig. n2] with the solid curve. The ratio of b$jb p indicates 
that the 4> meson interaction size, R™*, is smaller than 
that for the p meson: 



R l l 



Hint 
11 p 



= 0.87 ±0.08. 



(9) 
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* - CLAS Lukashin 

• - CORNELL Dixon 
■ - CORNELL Cassel 
▲ - HI Adloff 

O - SLAC Ballam, Q 2 = 
A - DESY Behrend, Q 2 = 
□ - BONN Besch, Q 2 = 
- CORNELL Berger, Q 2 = 
* - CLAS Anciant, Q 2 = 
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FIG. 12. The dependence of the t-slope, b^, on cAr. 
The electroproduction data CORNELL Dixon are from Ref. 
[17,18], CORNELL Cassel from Ref. [15], and HI Adloff from 
Ref. [34]. The photoproduction data BONN Besch are from 
Ref. [37], SLAC Ballam from Ref. [38], DESY Behrend from 
Ref. [42], and CLAS Anciant from Ref. [11]. 

A summary of the existing measurements of b^ to- 
gether with our results is shown in Fig. [ll]. Previous 
cj) electroproduction measurements are consistent with 
no Q 2 or cAr dependence JT^JT^] . We observe a low 
value of &0 « 2.2 GeV~ 2 , which, taken together with the 
values measured in photoproduction, shows a significant 
dependence on Q 2 . However, the Q 2 dependence of b^ 
can be explained by the implicit dependence of cAr on 



Q 2 (Eq. B7). This is shown in Fig.O where we plot the 
dep end ence of b^, on Q 2 using Eq. gand the relation in 



Eq. B7 at two values of W. The lower curve, at W = 2.3 
GeV, corresponds to our kinematics and connects our 
measurements with photoproduction values. The upper 
curve is closer to the Cornell kinematics. 
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FIG. 13. The dependence of the t-slope, b^, on Q 2 . Photo- 
production data BONN Besch are from Ref. [37], SLAC Bal- 
lam from Ref. [38], and CLAS Anciant from Ref. [11]. Elec- 
troproduction data CORNELL Dixon are from Ref. [17,18] 
and CORNELL Cassel from Ref. [15]. 

Because the value of cAr is smaller than the size of 
the nucleon, the scattering may be considered to be 
point-like. The application of QCD-inspired models, 
sensitive to the quark structure of the interacting me- 
son and nucleon, should provide an interesting inter- 
pretation of the observed 6(At) and b(Q 2 ) dependen- 
cies. It has been argued that with increasing Q 2 the 
radius of the virtual vector meson will shrink [jlj, and 
a corresponding decrease of b should be observed. At 
large enough Q 2 , quark models predict the de- 

crease of the tr ansverse d imension of the vector meson 
as r v - r h M/y/M 2 + Q 2 . The mass scale M represents 
a typical hadronic mass scale, which might be as small 
as the vector meson mass, e.g. = 1020 MeV, but 
is likely to be large compared to the Q 2 values of this 
experiment. Even though we do not need to invoke an 
explicit Q 2 dependence to describe our data, we note 
that the effects of transverse size and fluctuation times 
are not easily separated, especially when fine binning is 
prohibitive due to limited statistics. 
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IV. SUMMARY 



APPENDIX A: NOTATION 



The electroproduction of the 0(1020) vector meson 
was measured for Q 2 from 0.7 to 2.2 GeV 2 , W from 
2.0 to 2.6 GeV, and At from 1.8 to 4.0 GeV" 1 (cAr 
from 0.35 to 0.79 fin). A sample of 197 0(1020) mesons 
was accumulated for the exclusive reaction of ep — > e'p'tfi 
with the CLAS detector in Hall B at Jefferson Lab. 

(i) Taken together with the world data sample, we find 
a consistent picture of <j> production on the proton. Yet 
the scarcity of 4> data do not permit a precise quantita- 
tive description of the production process. 

(ii) We observe the expected decrease of the slope b$ 
of da/dt 1 on the formation length cAr below 2 fm. The 
rate of the decrease is similar to that in p meson pro- 
duction, but with a lower asymptotic value. Using a 
simple geometric model, the data show that the interac- 
tion size of 4> mesons with a proton is smaller than for p 
mesons. 

(iii) The <f> production cross section measurement adds 
new information at low values of Q 2 and W. The cross 
section dependence on Q 2 is qualitatively reproduced by 
pomeron exchange models. The cross section depen- 
dence on W as W°- 2±01 at Q 2 = 1.3 GeV 2 was de- 
termined by comparison to <fi production at HERA after 
correcting for threshold effects. This dependence is the 
same as observed in photoproduction. 

Additional electro- and photoproduction data from 
CLAS are currently being analyzed and will increase the 
overall qualitative and quantitative understanding of the 
physics that underlies vector meson production. 



We denote the four-momenta of the incident and scat- 
tered electron by p e and p e > , the virtual photon by q = 
Pe — Pe' , and the target and recoil proton by p p and p p > . 
Each four-vector can be written as (E, p) with appropri- 
ate subscripts. We use the common notation for Lorentz 
invariants: Q 2 = —q 2 > 0, v = q-p p /M p (M p is the mass 
of the proton) , the squared hadronic center-of-mass en- 
ergy W 2 = (q + Pp) 2 , and t = (p p — p P ') 2 is the four- 
momentum transfer to the target. The above-threshold 
momentum transfer is given by t' — t—t m i n (Q , W) < 0, 
where —t m i n is the minimum value of — t for fixed kine- 
matics. 




FIG. 14. Space-time picture of the yvP scattering through 
the conversion of the virtual photon into the virtual cj> meson 
inside the target proton. 



APPENDIX B: GEOMETRIC MODEL 
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We describe a qualitative picture of vector meson 
diffractive scattering within a simple geometric model. 
A sketch of the process is shown in Fig. |lj. The virtual 
photon is converted into the virtual vector meson (of ra- 
dius ry), which diffractively interacts with the proton (of 
radius r^) during a formation time At. Differential elas- 
tic cross sections are closely related to the charge form 
factors F(t) of colliding hadrons at high energy p9| , ^5[ . 
For small values of t, the form factors are related to the 
charge radii (r 2 ), via 



F(t) 



l + -{r 2 )t + 0(t 2 ) 



(Bl) 



For hadron-hadron elastic scattering (g9|, the cross sec- 
tions depend exponentially on t: 



da/dt 
(da/dt) t=0 



,M 



(B2) 
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Comparison of equations Bl and B2, and noting that 
the cross section is proportional to the square of the 
form factor, leads to a relationship between the radius 
of interaction, Ry *, and the i-slope parameter b: 

b=\{R^) 2 . (B3) 
The radius of interaction can be written as 

(R^r^{rl) + (r 2 v(Q 2 )), (B4) 

where r/j and ry are the radii of the nucleon and vector 
meson, respectively. 
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15. The i-slope parameter dependence on cAr for 
selected photo- and electroproduction data of p mesons. The 
data show a clear decrease of b with decreasing cAr below 2 



fm. The curve is a fit to Eq. B5 The photoproduction data 
SLAC Ballam are from Ref. [38] and Fermilab Francis from 
Ref. [46]. The electroproduction data CORNELL Cassel are 
from Ref. [15]. 

Because of the virtuality of the vector meson, the in- 
teraction region should also decrease if the formation dis- 
tance is less than the size of the nucleon (cAt ~ 2r^ rs 2 
fm) . A representative sample of the large body of p data 
shown in Fig. |l5| suggests the following phenomenological 
parameterization for the i-slope dependence on cAr: 



&(cAr) = 1 (1 - e - cAr ' 2rh ) (Ry 11 )' 



(B5) 



A two-parameter fit to Eq 
pendence of r„ on Q 2 , yielc 



B5, ignoring any explicit de- 



b P (cA T ) 



(7.86 ±0.26) 



1 



-cAr/2(0.78±0.05) 



(B6) 



However, Eq. Bl also has an indirect dependence on 
) 2 through cAr. At fixed W, we can write Eq. [I] as 



cAt 



c(W 2 



M 2 V 



Q 2 ) 



M p (Q 2 +M 2 ) 



(B7) 



|2 , using this 
for p data at 



Thus, we can plot Eq.[B5| as a function of 
expression for cAt. This is shown in Fig.[l 
the fixed value of W = 2.6 GeV |J|. Thus, we see that 
most, if not all, of the variation of the slope parameter b 
can be accounted for by changes in the fluctuation time. 

For the kinematics of this experiment, cAr (w 0.5 
fm) is small compared to the size of the nucleon, so we 
expect the fluctuation time factor to be significant for 
our d> data. 




O - SLAC Ballam, W » 2.6 GeV 
■ - CORNELL Cassel, W = 2.6 GeV 



-0.5 0.5 



1.5 



2.5 



3.5 



4.5 



Q 2 (GeV 2 ) 



FIG. 16. The f-slope parameter dependence on Q 2 for the 
photo- and electroproduction of p mesons at W = 2.6 GeV. 
The data show a clear decrease of b with increasing Q 2 . The 
The photoproduction data SLAC 
The electroproduction data COR- 
NELL Cassel are from Ref. [15]. 



curve is a fit to Eq. B5 
Ballam are from Ref. 



with x 2 /DF = 2.08. 
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